Relativistic quasiparticle time blocking approximation. II. Pygmy 
dipole resonance in neutron-rich nuclei 
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■ Abstract 



Theoretical studies of low-lying dipole strength in even-even spherical nuclei within the relativis- 
tic quasiparticle time blocking approximation (RQTBA) are presented. The RQTBA developed 
recently as an extension of the self-consistent relativistic quasiparticle random phase approximation 
(RQRPA) enables one to investigate effects of coupling of two-quasiparticle excitations to collective 
vibrations within a fully consistent calculation scheme based on covariant energy density functional 
theory. Dipole spectra of even-even ^^'^Sn - ^^'^Sn and ^'^Ni - ^'^Ni isotopes calculated within both 
RQRPA and RQTBA show two well separated collective structures: the higher-lying giant dipole 
resonance (GDR) and the lower-lying pygmy dipole resonance (PDR) which can be identified by a 
different behavior of the transition densities of states in these regions. 

PACS numbers: 21.10.-k, 21.60.-n, 24.10.Cn, 21.30.Fe, 21.60.Jz, 24.30.Gz 
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I. INTRODUCTION 



Since the last decade experimental and theoretical investigations of the pygmy dipole res- 
onance (PDR) remain a challenge of nuclear structure physics. Due to recent develop- 
ments of experimental techniques, in particular, the high resolution nuclear resonance flu- 
orescence, valuable information has been obtained about low-lying dipole spectra of stable 
nuclei Q, [sj. Pioneering Coulomb dissociation experiments 4, 5] allowed to extend inves- 
tigations of the low-lying dipole strength to nuclei with large neutron excess. 
Properties of the pygmy dipole resonance obtained in the various experimental and theoreti- 
cal studies can provide unique information about nuclear structure. It has been pointed out 
that the pygmy dipole strength is closely related to static nuclear properties, for instance, to 
the neutron skin thickness and the nuclear symmetry energy ^,0]. Besides the pure nuclear 
structure interest, an existence of a PDR just above the neutron threshold can have signifi- 
cant astrophysical importance as it may increase the neutron capture cross sections essential 
for r-process nucleosynthesis simulations (see, for instance, j?] and references therein). 
So far the best description of the experimentally observed low-lying dipole strength in sta- 
ble medium and heavy nuclei is achieved in the calculations within the quasiparticle-phonon 
model (QPM) [s^. The first 1~ state, the position, the total strength and the fine structure of 
the PDR are described very well in the wide model space including up to three-phonon con- 
figurations [l, 2|. However, because of lack of self-consistency and the presence of adjustable 
parameters, the application of this approach to nuclei with large neutron excess remains 
questionable. Another successful tool to investigate the low-energy dipole response are ap- 
proaches based on the covariant energy density functional, first of all, the fully self-consistent 
relativistic quasiparticle random phase approximation (RQRPA) see also Ref. [l^ for a 
recent review. The RQRPA supplemented with the coupling to low-lying vibrations within 
the relativistic quasiparticle time blocking approximation (RQTBA) HI] in a fully consistent 
way enables one to reproduce the fragmentation of the giant dipole resonance as well as of 
the PDR and to describe the dipole strength of the low-energy part of the spectrum. The 
coupling of the single-particle motion to vibrations is known as an important mechanism for 
the damping of particle-hole excitations. Ta king into account this mechanism is necessary to 



describe the widths of giant resonances 12 



Jj, Q, Q. Since the PDR is substantially 



a nuclear surface phenomenon, its coupling to surface vibrations may be significant. 
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The present paper is considered to be a contribution to the understanding of properties of 
the pygmy dipole resonance in neutron-rich nuclei. We investigate the dipole response of 
even-even spherical nuclei with large neutron excess within the RQTBA developed recently 



in Ref. 



ll\ paying specia 



the previous calculations 



attention to the lowest part of the dipole spectrum. Relying on 
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17l |. which demonstrate very good agreement with data for the 



experimentally known nuclei, we retain the same calculation scheme and do not introduce 
any additional parameters or approximations. 



II. RQTBA CALCULATIONS: AGREEMENT WITH DATA AND PREDICTIVE 
POWER 



In Ref. [ll| the relativistic quasiparticle time blocking approximation has been developed and 
applied for calculations of the dipole response of medium-mass open-shell nuclei. The physi- 
cal content of this approach is the quasiparticle-vibration coupling model based on covariant 
energy density functional theory (CEDFT) and the relativistic QRPA (RQRPA) f^. To 
describe excitations in nuclei with even particle number the Bethe-Salpeter equation (BSE) 
containing both static and energy- dependent residual interactions has been formulated for 
the nuclear response function R{u) in the doubled quasiparticle space: 



where the indices ki run over the single-particle quantum numbers including states in the 
Dirac sea and indices r], t]', rj" numerate forward (+) and backward (-) components in the 
doubled quasiparticle space. The quantity w-^\(jj) describes the free propagation of two 
quasiparticles in the mean field between their interaction with the amplitude 

^kikiMk^i^) = K'X.fefca + i^lik^Mki^^) ~ ^Iiki,k2k3(^))^m'- (2) 

The static part of the^ingle-quasiparticle self-energy, determined by the CEDFT with the 
parameter set NL3 [18], is based on a one- meson exchange interaction with a non-linear 
self-coupling between the mesons. Pairing correlations, considered to be a non-relativistic 
effect and treated in terms of Bogoliubov's quasiparticles within the BCS approximation, 
are introduced into the relativistic energy functional as an independent phenomenologically 
parameterized term. The energy-dependent part of the single-quasiparticle self-energy is 



modelled by the quasiparticle-phonons coupling (QPC). The energy-dependent residual in- 
teraction is derived from the energy-dependent self-energy by the consistency condition 
and calculated within the quasiparticle time blocking approximation. This approximation 
means that, due to the time projection in the integral part of the BSE, the two-body propa- 
gation through states which have a more complicated structure than 2q ® phonon has been 



blocked 15|. In order to avoid double counting of the QPC effects a proper subtraction of the 
static QPC contribution from the residual interaction has been performed in Eq. ([2]), since 
the parameters of the CEDFT have been adjusted to experiment and include already essen- 
tial ground state correlations. In order to take the QPC into account in a consistent way, 
we have first calculated the amplitudes of this coupling within the self-consistent RQRPA 
with the static residual interaction V. Then, the BSE ([1]) for the response function has been 
solved in both Dirac-Hartree-BCS and momentum-channel representations. 
To describe the observed spectrum of an excited nucleus in a weak external field P as, for 
instance, an electromagnetic field, the microscopic strength function S{E) is computed: 

^(^) = -^ii- ^ Y.Pr.^MLMH(E + ^^)P^,^. (3) 

In the calculations a finite imaginary part A of the energy variable is introduced for con- 
venience in order to obtain a smoothed envelope of the spectrum. This parameter has the 
meaning of an additional artificial width for each excitation. This width emulates effec- 
tively contributions from configurations which are not taken into account explicitly in our 
approach. The dipole photoabsorption cross section aEi{E) is expressed via the strength 
function Sei{E) according to the well known formula: 

167r^e^ 

aEiiE) = ^^ESEiiE). (4) 



In Ref. Ill] numerical results for dipole spectra in tin isotopes and N=50 isotones have been 
presented, which were obtained within the approach described above and compared with 
RQRPA calculations and the available experimental data. The mean energies (ME) and 
the energy-weighted sum rule (EWSR) of the giant dipole resonance (GDR) in all nuclei 
under investigation agree very well with the experimental data, and this is provided by the 
self-consistent RQRPA. The coupling to phonons is responsible for the fragmentation of the 
resonance. In our fully consistent approach it does almost not affect the nice agreement 
achieved in RQRPA for the ME and the EWSR in the GDR region. 



Notice, that in stable nuclei not only the integral characteristics like the mean energy, the 
energy-weighted sum rule and the width, but also the general shape of the GDR are obtained 
in the RQTBA rather close to those observed in experiments demonstrating the new quality 



of description as compared to previous microscopic studies. Thus, from Ref. where 
an excellent agreement has been achieved with the data for the GDR in ^^^Sn, ^^''Sn, ^^Sr, 
^°Zr, we conclude that the main mechanisms of the GDR formation in medium mass and 
heavy nuclei are taken into account correctly and consistently in the RQTBA. Moreover, 
it turned out to be possible to obtain a very good description of the low-lying part of the 
dipole strength below the nucleon separation energy in the same calculation scheme. Strictly 
speaking, a correct comparison of the integrated pygmy dipole strength with experimental 
data in this region is still an open problem since the existing data, which result in most 
cases from the (7,7') measurements, are usually restricted from the top by the nucleon 
separation energies because of loss of experimental sensitivity. In theoretical calculations 
beyond (Q)RPA the pygmy mode is in most cases strongly coupled to the surface vibrations 
which causes the spreading of the low-lying strength with the isoscalar underlying structure 
in a broad energy region, for instance, 3-10 MeV for tin isotopes with medium neutron 
excess. So, at present the comparison with data is only possible for the lowest part of 
the pygmy strength measured below 8 MeV. For QRPA calculations the situation is even 
more difficult: in particular, in self-consistent RQRPA calculations only one or very few 
pronounced peaks, corresponding to the soft dipole mode, occur and accumulate the whole 
strength with the isoscalar underlying structure, so that they often do not overlap at all with 
the experimentally investigated region. However, the situation is different for the RQTBA 
approach. Within this model ll|| the integrated dipole strength below 8 MeV is in very 
good agreement with the data which are available for ^^Sr 19| and ^^^Sn Note, that 
the coupling to phonons is the essential mechanism which brings the strength to this region 
while the RQRPA yields no dipole strength at all in this energy range. 
Thus, the RQTBA has been successfully tested and demonstrated as a reliable approach for 
the description of excited states in stable nuclei in the GDR as well as in the PDR regions. 
Since it contains no adjustable parameters and uses only the universal relativistic mean field 
(RMF) parameters, the approach should also give reliable results for nuclei further away 
from the valley of stability. This expectation has been indeed confirmed by a comparison 
with recent pioneering experiments on Coulomb dissociation of i30,i32gj^ rpj^^ dipole 
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FIG. 1: Spectra of dipole excitations in ^^'^Sn, ^'^^Sn. The experimental data [J| are shown by the 
black circles with the error bars, the thin curves represent our RQTBA calculations with 1 MeV 
smearing, and the thick curves give the RQTBA strength distributions after the convolution with 
the detector response 111]. 

strength calculated within the RQTBA for the doubly magic ^^^Sn has been presented in 
Ref. 171], and its low-lying part has been analyzed in Ref. 2Q[. Calculations for ^^°Sn have 
been done and presented in Ref. However, a direct comparison of these results to the 
data of Ref. {4] is not possible, but rather requires a convolution of the calculated strength 



with the detector response matrix 



2l[ |. For this purpose we have first recomputed the dipole 



strength in these two nuclei with the smearing parameter 1 MeV which corresponds to the 
experimental energy bin size (note, that in Refs. [ll, 17, 2o| the calculations were made 
with 200 keV smearing). The results are presented by the thin solid curves in Fig. [H The 
thick solid curves give the results obtained after the convolution with the detector response. 
The black circles with error bars are the experimental data. After the convolution, the 
RQTBA dipole strength distributions reproduce the data very well, although the major 
part of the RQTBA pygmy strength resides below the neutron threshold, where the dipole 
strength is not detectable in the experimental set-up of Ref. Moreover, it turns out that. 
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after the convolution, the data can be similarly well reproduced by quite distinct initial 
dipole strength distributions; e.g. by those obtained within QRPA based on the Skyrme 
functional [6]. Thus, one may conclude that the data are only weakly sensitive to the details 
of the original strength distribution and, for a better understanding of the PDR phenomenon 
in these nuclei more precise data are needed. 



III. PYGMY DIPOLE RESONANCE: GENERAL FEATURES 

Ref. 



201 ] reports about detailed investigations of low-lying dipole strengths in nuclei with 
closed shells within the relativistic time blocking approximation (RTBA). It has been con- 
firmed that the PDR arises from a resonant oscillation of the neutron skin against the isospin 
saturated proton-neutron core, with the corresponding relativistic RPA state characterized 
by a coherent superposition of many neutron particle-hole configurations. Moreover, this pic- 
ture remains essentially unchanged when the particle-vibration coupling is included within 
the RTBA. The effect of two-phonon admixtures in nuclei with closed shells results in a weak 
fragmentation and a small shift of the PDR states to lower excitation energies, while the 
one-phonon character and the underlying structure of the states are not noticeably modified 
by the coupling to low-lying surface vibrations. 
However, calculations performed in Ref. 



11| have shown that the presence of pairing cor- 
relations at least in one of the nucleonic subsystems may change this conclusion crucially. 
First of all, these correlations lead to a diffuseness of the Fermi surface and, thus, increase 
the number of possible configurations of 2q ® phonon type. But the main reason for the 
stronger coupling effect is the considerable lowering of the energies and/or the increased 
transition probabilities of the lowest 2'^ states due to pairing correlations of the superfluid 
type. In spherical open-shell medium-mass nuclei highly collective first 2+ states appear at 
energies around 1 MeV (and they are usually well reproduced in RQRPA) whereas in magic 
nuclei and often in nuclei near shell closures they appear much higher, at about 3-4 MeV and 
usually have a considerably reduced strengths. This has important physical consequences 
resulting in strong configuration mixing if low-lying vibrational states are present and the 
fragmentation caused by this mixing is the stronger the lower the excitation energies and 



the higher the transition probabilities. This is confirmed by the calculations of the Ref. 



In the RQTBA calculations we can, in particular, analyze admixtures to the pygmy mode 
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which are associated with the collective vibrational motion of the neutron or proton excess 
against the isospin saturated core. This collective state, whose collectivity is directly depen- 
dent on the neutron excess (i.e. the |N-Z| value), appears in the RQRPA as a pure mode with 
a pronounced dominance of the neutron or proton component of the transition densities at 
the nuclear surface, which, in turn, is subject to vibrational motion. As a consequence, the 
pygmy mode is not observed as a pure mode, but always couples to the surface vibrations. 
The results obtained in Ref. llj for Z=50 and N=50 nuclei with medium neutron excess 



also show that coupling to low-lying vibrational states is the main mechanism which brings 
the dipole strength to the low-energy region far below the neutron threshold, whereas in 
pure RQRPA only one collective peak (or several peaks in nuclei with larger neutron excess) 
around the threshold occurs. This peak, which arises as a result of a pushing-down effect 
and the coherence of the unperturbed two-quasiparticle states by the isoscalar part of the 
effective meson-exchange interaction, then spreads over a broad energy range due to phonon 
coupling. As a result, one finds its fragments well below and above its original position in 
close resemblance of the spreading mechanism of the GDR. 



IV. LOW-LYING DIPOLE STRENGTH IN NEUTRON-RICH TIN ISOTOPES 



In Ref. 



ll| the dipole response in the tin isotopic chain has been investigated up to ^^°Sn 



(with neutron number N=80), while results for the doubly magic nucleus ^^^Sn are given in 
Refs. p/7|,l20|. Here we extend these studies of the dipole response in even-even tin isotopes to 
calculation of the very neutron-rich isotopes ^^^Sn, ^^^Sn, ^^^Sn, ^^'^Sn. The results are shown 
in Figs. [2] and [31 including also results obtained previously for ^^''Sn and ^^^Sn for comparison. 
The right panels of the figures display the photo absorption cross sections, calculated with 
the usual isovector dipole operator, using a 200 keV smearing. The left panels enlarges the 
low-lying parts of the corresponding spectra, calculated with a smaller imaginary part for 
the energy variable (20 keV), in order to pronounce the fine structure of the spectra and to 
even identify individual levels in this energy region. Calculations within the RQRPA and 
RQTBA are shown by the dashed and solid curves, respectively. The arrows denote the 
positions of the experimental and the theoretical neutron thresholds, showing reasonably 
good agreement. 

Fig. [2] represents the case of the doubly magic nucleus ^^^Sn and its neighboring even-even 
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FIG. 2: The calculated dipole spectra for I30.i32,i34gj^ isotopes. Right panels: photoabsorption 
cross sections computed with the artificial width 200 keV. Left panels: the low-lying portions of 
the corresponding spectra in terms of the strength function, calculated with 20 keV smearing. 
Calculations within the RQRPA are shown by the dashed curves, and the RQTBA - by the solid 
curves. Arrows denote the experimental and the theoretical neutron separation energies. 



nuclei ^^°Sn and ^^^Sn. Comparing these results to those obtained in Ref. ll| for ^^^Sn 
and ^^°Sn and with the more neutron-rich isotopes shown in Fig. [3l one observes that the 
fragmentation of the FDR is much weaker in nuclei around shell closure. This confirms our 
general discussion about fragmentation above since in the nuclei close to the magic numbers 
the excitation energies of the low-lying 2^ states are noticeably higher and the transition 
probabilities are reduced compared to the situation in open-shell nuclei. This means that 
the respective matrix elements of the phonon vertices have also smaller values, leading to a 
weaker nucleon-vibration coupling. 
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FIG. 3: Same as in Fig. El but for ise.iss.uogj^ isotopes. 

We further observe that the widths of the GDR are also reduced in magic nuclei as compared 
to nuclei with opened shells. While typical width values are 5-5.5 MeV [11] in agreement 
with data and R(Q)TBA calculations, they are reduced to 3.4 and 4.0 MeV in the R(Q)TBA 



calculations for the doubly magic ^°°Sn ll| and ^^^Sn isotopes. This reduction has, obviously, 
a similar origin as the reduction of the PDR fragmentation discussed above, confirming a 
strong resemblance between the GDR and PDR features. 

The smooth decrease of the GDR mean energy with the nuclear mass number A, predicted 
by the well known systematic rule {E) ^ 80A^^^^, is well reproduced in our RQTBA calcula- 
tions. Our results for the GDR and the PDR mean energies are also in a good agreement with 
the results of the RQRPA studies of Ref. [22| for tin and nickel isotopes obtained with the 



density-dependent meson-exchange forces (DD-MEl). Visually, the whole dipole spectrum 
including both the GDR and the PDR regions shifts towards lower energies with increasing 
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mass number. However by closer inspection, we observe that by increasing only the neu- 
tron number and, thus the isospin asymmetry, the fragmented pygmy resonance moves down 
quicker than the GDR. This is illustrated in Fig. HI Whereas in ^^^Sn and in the isotopes 
with even smaller neutron excess the PDR is barely visible in the cross sections, for more 
neutron-rich nuclei, starting from ^^''Sn it appears as a pronounced structure which is well 
separated from the GDR, in both the RQRPA and RQTBA. This can be understood in 
terms of the well known Brown-Bolsterli model {23! where an increase of collectivity and 
also of the strength of the lowest solution is directly connected to its shift towards lower 
energies. This qualitative observation about the separation of the PDR from the GDR in 
neutron-rich nuclei will be confirmed by the analysis of transition densities in the energy 
region around 10 MeV which we present below (see Figs. O O [71 [8]) . In Fig. [3] one can see 
how the GDR and PDR strengths develop with a further increase of the neutron excess. The 
PDR continues its separation moving to lower energies, and its fragmentation also becomes 
stronger as well as the fragmentation of the GDR. The latter is related to the considerable 
lowering and increasing collectivity of the first 2+ states. For example, the RQRPA predicts 
the first 2+ state at 0.6 MeV in ^^''Sn indicating a trend to a static quadrupole deformation 
in this nucleus. Indeed, the systematic deformed relativistic Hartree-BCS calculations with 
the NL3 force of Ref. 2J] indicate a phase-transition from spherical to axially deformed 
shapes in this region. 

Before we study the structure of the PDR in neutron-rich nuclei in some details, we would 
like to make few comments about their astrophysical importance. As has been pointed out by 
Goriely and co-workers [?], an enhancement of the dipole strength at energies just above the 
neutron thresholds can increase the neutron capture cross sections relevant for simulations 
of the astrophysical r-process. From the tin isotopes studied here, the double-magic ^^^Sn 
and the even more neutron-rich nuclei are close to the r-process path. As can be seen in 
Fig. [21 our calculation predicts a sizable contribution of the PDR in ^^^Sn just above the 
neutron threshold which is, however, relatively high (about 7 MeV) in this double-magic 
nucleus. The presence of such PDR strength increases the neutron capture cross section 
on ^^^Sn compared to the usual treatment which models the low-lying dipole strength by 
a Lorentzian fit to the GDR strength [25|]. As explained above the location of the PDR 
is pushed down in energy with increasing neutron excess. However, the neutron threshold 
energy decreases even faster in the tin isotopes beyond ^^^Sn. As a consequence the presence 
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FIG. 4: The mean energies for the GDR (filled symbols) and the PDR (opened symbols) in tin 
isotopes calculated within the RQTBA. 

of the low-lying PDR strength, predicted by our RQTBA calculations, is not expected to 
increase the corresponding neutron capture cross sections as it is located noticeably above 
the respective neutron thresholds. To the contrary, the lack of dipole strength just above the 
neutron thresholds, as predicted by our RQTBA calculations, might result in lower capture 
cross sections than obtained from the parameterized Lorentzian strength distributions. An 
answer to this question will be given by statistical model calculations which will compare 
cross sections obtained from our microscopically derived dipole strengths with those using 
the standard Lorentzian approximation. These studies are in progress. 
The rather strong first dipole states found in our RQTBA calculations for the nuclei 
I36,i38,i40gj-^ can be interpreted as the prototypes of the 2~^^3~ states which are usually 
observed in dipole spectra and always appear, for instance, in QPM calculations (see, for 
instance, Ref. [l|). Since in the present version of the RQTBA we confine ourselves to 
2q phonon configurations, the position of each peak, if it would be unperturbed by the 
static residual interaction, is determined by the sum of the two quasiparticle energies and 
the phonon energy, whereas in the QPM it is determined by the sum of two phonon energies 
which is, as a rule, considerably lower in energy and has correspondingly often an enhanced 
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transition probability. A similar effect should take place in the two-phonon version of our 
model where instead of two unperturbed quasiparticles their bound state, a phonon, is used 
as an elementary excitation (see 15| for the details). In the present case the 2q ® phonon 
configurations in nuclei with medium neutron excess, e.g. configurations consisting of two 
quasiparticles with total angular momentum and parity 3^ and a 2+ phonon, have higher 
energies and reduced probabilities, so that they are not noticeable in the lowest part of the 
spectrum. In more neutron-rich, nuclei, however, the 2+ states are pushed down to such low 
energies and relatively strong transition probabilities that their matrix elements are large 
enough to form a low-lying 2q ® phonon configuration with the sizable strength. We note 
that the same explanation is applicable to the case of neutron-rich nickel isotopes to be 
discussed below (see Figs. El |9]). 

For a detailed inspection of the structure of prominent peaks in the RQRPA and RQTBA 
dipole strength distributions we have calculated their respective transition densities. For 
comparison we focus on two nuclei: the stable isotope ^^^Sn with a rather small neutron excess 
(Figs.[5l[6]) and ^^°Sn with a large isospin asymmetry (Figs. [71 [8]). In the top panels of figures[5] 
and[7]we display the low-lying dipole strength distribution, calculated within the RQRPA and 
RQTBA. The middle and bottom panels of these figures exhibit the neutron (solid curves) 
and the proton (dashed curves) RQRPA transition densities for the most pronounced peak 
of the PDR region (middle) and the lowest peak of the GDR region (bottom) as obtained 
within the RQRPA approach. The choice of the latter peak indicates approximately the 
energy region where the PDR can be separated from the low-energy tail of the GDR. 
Fig. [5] shows that in the RQRPA model nearly the entire pygmy strength in ^^^Sn is concen- 
trated in a single very strong peak localized at 8.94 MeV. The neutron component obviously 
dominates the transition density at the surface, however, with a noticeable proton admix- 
ture (middle panel of Fig. [5l). In contrast, the corresponding panel of Fig. [7] implies that the 
proton component is fully suppressed at the nuclear surface in the pygmy mode of ^'^''Sn, 
corresponding to the state at 7.18 MeV. The bottom panels of Figs. [5] and [7] show that the 
structure of the RQRPA peaks changes drastically when the excitation energy increases and 
the lower part of the GDR region is reached. Already the peaks at 11.78 MeV in ^^^Sn and at 
10.94 MeV in ^''^Sn exhibit a completely different relative behavior of the neutron and proton 
components. The two components are still in phase in the nuclear interior for ^^^Sn, while 
they are not for ^^°Sn. For ^^^Sn neutrons and protons are out of phase at the surface with 
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FIG. 5: The low-lying part of the dipole spectrum in ^^^Sn, calculated in RQRPA (dashed curve) 
and RQTBA (solid curve) models (upper panel). RQRPA transitional densities for the most inten- 
sive peaks at 8.94 MeV (middle) and 11.78 MeV (bottom). The dashed curves represent the proton, 
and the solid curves - the neutron transitional densities. See text for the detailed explanation. 

protons even dominating, whereas in ^^°Sn protons do not contribute to the surface motion. 
Thus the bottom panels indicate that the states at these moderate excitation energies are no 
longer of pygmy type, but do also not yet exhibit the typical GDR structure with protons 
and neutrons oscillating against each other. 

In Figs. [6] and [8] we display the dipole transition densities for selected states in ^^^Sn and 
^^'^Sn, as calculated within the RQTBA. These states arise from the coupling of the RQRPA 
states to 2^ 3~, 4"'", 5^, 6^ phonons with energies below 10 MeV. As it has been shown in 



Refs. 171, l20|] for the case without pairing correlations, the RTBA transition densities obey 
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FIG. 6: The RQTBA transition densities for the low-lying dipole states in ^^^Sn. The dashed 
curves represent the proton, and the solid curves - the neutron transition densities. 

more general normalization condition than the pure RPA states. This normalization condi- 
tion is generalized straightforwardly to systems with pairing correlations of the superfiuid 
type. Formally, this condition in the RQTBA has the same form as in the RTBA, but the 
respective summation is performed over the doubled quasiparticle space. In the notations 



of Ref. the deviation of the quantity J2 ^l'^JI(fcifc2)l^ (where '^^^(^kik2) matrix 

(kik2)r) 

elements of the transition densities of the state /i) from unity indicates quantitatively the 
effect of the coupling to phonons of the state /i. Notice, that the amplitudes of oscillations 
are reduced in the RQTBA as compared to the RQRPA as reflected in Figs. [6] and [HI These 
figures give a detailed look into the underlying structure of the RQTBA dipole spectrum. In 
particular, we observe that in the case of strong configuration mixing, like in the RQTBA, 
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FIG. 7: Same as in Fig. El but for ^^OSn. 

both nuclei develop a transitional region with states which have neither pure isoscalar nor 
isovector underlying structures. This is in strong contrast to the RQRPA which does not 
show such a smooth transition range. For example, except for two relatively weak excita- 
tions around 9 MeV, the RQRPA spectrum for ^^°Sn does not contain states in the energy 
range between 8 and 11 MeV. Thus we conclude that the RQRPA predicts a clear separation 
between the PDR and GDR regions. This separation is, however, overcome in the RQTBA 
due to strong mode mixing giving rize to the formation of a transitional region between PDR 
and GDR containing many weak states with mixed structure. 

Table [T] presents distributions of two-quasiparticle configurations which contribute to the 
selected low- lying dipole states in ^^^Sn: in the first column we display the configurations 
forming the RQRPA state at 8.94 MeV which accumulates almost the entire strength with 
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FIG. 8: Same as in Fig. El but for ^^^Sn. 

isoscalar underlying structure within the RQRPA. From Table H] one can see that this is a 
collective state since many 2q configurations contribute with comparable strengths. From 
this distribution one can also conclude that in this nucleus, although neutron configurations 
dominate the pygmy mode, admixtures of proton configurations are still sizable, since the 
neutron excess is not yet large. As it has been discussed above, in the RQTBA this state is 
strongly coupled to phonons and fragments within the energy region between 4 and 10 MeV 
into many states with isoscalar underlying structure. The second, third and fourth columns 
of Table [T] represent the 2q contributions to the states at 6.44 MeV, 9.20 MeV and 10.23 MeV, 
respectively. Note, that the total percentage of the 2q configurations is very much reduced 
with respect to the 100 % normalization which one has in the RQRPA. This reduction from 
unity is a measure of the coupling strength to phonons for a certain state. We find that for 
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TABLE I: Distribution of two-quasipaxticle configurations contributing to the low-lying dipole 
states in ^^^Sn calculated within the RQRPA and the RQTBA. 
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the considered states in this region deviations are about 50-70 % that indicates very strong 
coupling. 

V. RQTBA STUDIES FOR NEUTRON-RICH NICKEL ISOTOPES 

In this section we present the results of our studies of the dipole strength in neutron-rich 
nickel isotopes. One of our motivations for choosing these isotopes is the recent experi- 
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mental work [26] where the low-lying dipole strength in ^^Ni has been studied by Coulomb 
dissociation, observing a pronounced and broad bump in the excitation spectrum with a 
centroid between 10 and 11 MeV. Our second motivation stems from nuclear astrophysics as 
the double-magic nucleus ''^Ni is considered the essential waiting point at the = 50 shell 



closure for certain r-process nucleosynthesis simulations [25|. 

Analogously to the presentations of the tin isotopes in Figs. [2] and [3], we display the RQRPA 
and RQTBA dipole spectra for the even nickel isotopes ^^Ni - '''^Ni in Figs. M and [TDl The 
overall pattern of the dipole spectra (right panels), including the fragmentation of the GDR 
and of the low-energy part of the spectra, is quite similar for ^^Ni and ™Ni, which differ only 
by two neutrons in the l(?9/2 orbit within the independent particle model. The similarity 
of the fragmentation effects is caused by the fact that the first 2+ states, which are mainly 
responsible for the fragmentation, have comparable energies and transition probabilities in 
the RQRPA for both nuclei. In '^2,74,76]\jj ^j^g PDR and GDR are stronger fragmented due to 
lowering of the energies and an accompanying increase of the transition probabilities of the 
first 2+ states. As expected, in the double-magic nucleus ^^Ni the fragmentation is reduced 
again. 

Let us discuss our results for ^^Ni in more details. For the low-lying part of the dipole 
spectrum (Fig. [9](a)), the RQRPA calculation (dashed curve) predicts one rather pronounced 
peak at 9.52 MeV, accompanied by three relatively weak states. The transition density 
of the main peak has the isoscalar underlying structure typical for the pygmy mode (see 
Fig- [m below). The enlarged picture of the low-energy part of the spectrum, calculated with 
the smaller smearing parameter (upper panel of Fig. [TT!) indicates that the next RQRPA 
structures appear only at around 14 MeV. The transition densities of these states do not 
have an isoscalar underlying structure in the nuclear interior, but behave already as the 
low-energy tail of the GDR. Thus, in this nucleus, as well as for the heavier nickel isotopes 
(Figs. [9] and [TOl) . we observe the pronounced separation in the RQRPA dipole strength 
distribution between the PDR and the GDR regions, as already familiar from the neutron- 
rich tin isotopes. 

When coupled to phonons the pygmy mode, which is localized in the RQRPA dipole spec- 
trum for ^^Ni between 9 and 10.5 MeV, spreads over the broad energy region from 7.5 to 
about 12 MeV. The situation is very similar in the other investigated nickel isotopes. The 
corresponding RQTBA strength distributions in Figs. [9], [10] show many states with compa- 



19 




FIG. 9: Same as in Fig. [21 but for 68,70,72]v^j isotopes. 

rable strength in the energy region below 12 MeV. One of the lowest sharp peaks in ®^Ni - 
''^Ni is supposed to have a 3" ® 2+ 2g ® phonon nature. Hence the energy of such a state 
is expected to be pushed down in the two-phonon version of the RQTBA. Among the single 
peaks at low energies in ^^Ni one can see a concentration of strength centered at 10.60 MeV 



26j- In Fig. [Hthe 



which might correspond to the experimentally observed peak in Ref. 
transition densities of selected RQTBA states in the low-lying region are displayed. One can 
see that, in the case of ^^Ni only the few lowest states exhibit a pure isoscalar underlying 
structure. In the energy region between 10.5 and 13.5 MeV, for which the RQRPA does not 
produce any dipole strength, the RQTBA predicts a large amount of states with transitional 
behavior: their density distributions have neither pure isoscalar nor isovector underlying 
structure, however, the isoscalar character dominates below 12 MeV. 

Astrophysically, ''^Ni plays an important role for the matter flow beyond the = 50 shell 
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FIG. 10: Same as in Fig. El but for ^4,76,78^^^ isotopes. 

closure in certain r-rocess scenarios. This has been demonstrated in Ref. reporting 
the first measurement of the ^^Ni half hfe which turned out to be noticeably shorter than 
predicted by the global models and had a decisive influence of the local r-process abundance 
distribution around the = 50 nuclides. Our focus here is the question whether our 
microscopic approach predicts an enhanced dipole strength at energies just above the neutron 
threshold which might affect the neutron capture cross section to ^^Ni. As can be seen in 
Fig. [To], our RQTBA calculation indeed predicts the existence of a strong pygmy strength 
in ''^Ni, which, however, resides more than 2 MeV above the neutron threshold and hence is 
expected to have only a mild effect on the capture cross sections. Nevertheless, statistical 
model calculations employing our RQTBA dipole strength distributions are in progress for 
all nickel isotopes. 

As a more general remark we note that our RQTBA results should be reliable enough to 
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FIG. 11: Same as in Fig. [5l but for the ^^Ni. 

check whether the Lorentzian approximation to the dipole strength function as imphed from 
data for stable nuclei and often employed to derive the relevant neutron capture cross sections 
for r-process nuclei is also valid for nuclei with large neutron excess. Such studies are also 
in progress. 

VI. SYSTEMATIC OF THE PYGMY DIPOLE STRENGTH 

This section is devoted to a study of the systematics of the total PDR strength as a function 
of neutron excess. Obviously, the crucial point of such an investigation is the determination 
of the energy at which the PDR and GDR regions separate; i.e. the energy at which the 
underlying structure of the states in the dipole strength function change from an isoscalar 
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FIG. 12: Same as in Fig. El but for the ^^Ni. 

to an isovector character. We have determined this energy by a careful inspection of the 
transition densities of the individual states and their underlying structure. Having deter- 
mined in this way the upper energy limit, we have calculated the total PDR strength by 
integration of the strength function obtained within the RQRPA and the RQTBA with a 
smearing parameter of 20 keV. The left panel of Fig. [13] displays the integrated RQRPA and 
RQTBA pygmy strength ^ B{E1) f for the tin isotopes versus their neutron number. In 
the right panel one can see the integrated pygmy strength calculated within the RQTBA for 
tin and nickel isotopes as well as for ^osp]^ functions of the squared asymmetry parameter 
a={N - Z)/A. 

For the tin isotopes the total pygmy strength generally increases with neutron number, both 
in the RQRPA and RQTBA. One can see from Fig. [T3l(a) that the difference between the 
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FIG. 13: Integrated dipole strength with isoscalar underlying structure below GDR (pygmy 
strength) calculated within the RQTBA. Left panel (a): pygmy strength in tin isotopes, calcu- 
lated within the RQRPA (open circles) and the RQTBA (solid squares), versus the neutron number. 
Right panel (b): pygmy strength in nuclei versus the squared asymmetry parameter a = [N — Z)/A 
in Sn isotopes (solid squares), Ni isotopes (open squares), and ^'^'^Pb (triangle). 

RQRPA and the RQTBA results for the integrated pygmy strength is not sizable. This 
confirms our general suggestion that the phonon coupling causes fragmentation and redistri- 
bution of the pygmy mode without breaking its integrity, as it happens with the GDR. Some 
difference between the RQRPA and the RQTBA total pygmy strength may be explained by 
uncertainties in determining of the upper integration energy limit. Being well determined 
in the RQRPA, this quantity is only approximate in the RQTBA case because for states 
with the transitional underlying structure we rely only on our visual impression deciding if 
a certain state has more isoscalar or more isovector underlying structure. Nevertheless, in 
the present study we have determined the value 10 MeV as the upper integration limit for 
the pygmy strength in tin isotopes, 12 MeV - for nickel isotopes and 8 MeV - for ^''^Pb. 
However, the behavior of the integrated pygmy strength shown in Fig. [13] is not purely linear: 
our RQTBA calculations predict a nearly constant pygmy strength for the tin isotopes 
while the RQRPA pygmy strength shows also a deviation from the linear 
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behavior within this isotopic range. Note, that these isotopes have in common that the 
the last (partially) occupied neutron subshell is the l/iii/2 orbit, which is closed at N=82. 
In the right panel, where we present the RQTBA pygmy strength versus the asymmetry 
parameter, one can see that the total pygmy strength in the neutron-rich nickel isotopes 
68]\]'j_78]\j-j (demonstrates a very similar behavior: it remains almost constant for these isotopes. 
Interestingly, for these nuclei the neutron lgg/2 subshell is also being filled up closing the 
N=50 shell. In addition, the lgg/2 orbit and the l/iii/2 orbit are intruder orbits in the Ni 
isotopes and Sn isotopes, respectively. 

Another interesting observation is that the absolute values of the pygmy strength in tin and 
nickel isotopes with the same value of the asymmetry parameter are roughly close to each 
other and even to the value deduced for ^osp]^ from the calculations of Ref. 17|]. We observe 
quite a big difference between the total pygmy strength in Sn and Ni chains only at > 0.06 
that can be explained by approaching the shell closure N = 50 in the nickel chain. Certainly, 
to make a general statement about the dependence of the pygmy strength on the asymmetry 
parameter, one should perform a more wide investigation of different isotopic chains, but 
our observation agrees, in general, with the observation of Ref. ^ based on the systematic 
of the available experimental data. 

Notice, that in the present studies we use the NL3 interaction which is known to overestimate 
systematically the neutron skin thickness {28]. The reason is that in the NL3 parameter set 
the meson coupling in the isovector channel has no density dependence and is parameterized 
by only one coupling constant. As a result, the total pygmy strength in nuclei with large 
neutron excess can be respectively overestimated as compared, for instance, to the values 
obtained in the RQRPA calculations 221] with DD-ME2 interaction presented in Ref. 29|. 
From our systematic studies we conclude that in the models based on the CEDFT the 



total pygmy strength can deviate from the 
the asymmetry parameter. In Refs. [s^. 



inear dependence on the neutron number or on 



3l| it has been found within a non-relativistic 



model with Landau-Migdal interaction, that the coupling to collective vibrations produces a 
noticeable deviation of the pygmy dipole strength from the linear dependence on the neutron 
number in calcium isotopes, in agreement with the experimental data of Ref. jsO]. From our 
present investigation for the tin chain we conclude that this deviation can take place around 
the shell closures with intruder orbits both for the RQRPA and the RQTBA systematics. 
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VII. SUMMARY 

The relativistic quasiparticle time blocking approximation developed recently is applied 
to a systematic study of dipole spectra in even-even neutron-rich tin and nickel nuclei with 
a special focus on the low-lying part of the dipole strength. 

First of all, a very reasonable agreement of the RQTBA results with recent measurements for 
^^°Sn, ^^^Sn jj], which determined the dipole strength above the neutron thresholds, is found. 
The RQTBA is also applied to calculations for nuclei with large isospin asymmetry paying 
special attention to the fragmentation of the GDR and the PDR modes due to coupling to 
phonons. It has been found that the fragmentation is considerably reduced in doubly-magic 
nuclei and in some cases in neighboring even-even ones. The connection of this effect with 
the properties of the lowest 2"*" states has been discussed. 

The transition densities of the low-lying dipole states have been computed within the RQRPA 
and the RQTBA. It has been shown that the pygmy mode with the isoscalar underlying 
structure spreads in the broad energy region over many states whose underlying structure 
demonstrates a similar type of behavior which helps to identify the fragments of the pygmy 
mode and to separate them from the low-energy GDR tail. With help of this analysis the 
integral pygmy strengths in the investigated nuclei have been computed. It has been shown 
that the calculated integral pygmy strength demonstrates a nearly linear behavior versus the 
asymmetry parameter with some deviations from this behavior near the shell closures. 
As a next step to astrophysical applications, the RQTBA dipole strength distributions will 
be used in statistical model calculations to derive neutron capture cross sections for r-process 
simulations. 
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